Evaluating transport coefficients in real time thermal field theory 
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Transport coefficients in a hadronic gas have been caicuiated eariier in the imaginary time for- 
muiation of thermai fieid theory. The steps invoived are to reiate the defining retarded correlation 
function to the corresponding time-ordered one and to evaluate the latter in the conventional per- 
turbation expansion. Here we carry out both the steps in the real time formulation. 



o. 

(N . I. INTRODUCTION 

Thermal quantum field theory has been formulated in the imaginary as well as real time [1-5] . For time independent 
| /" ( 1 quantities such as the partition function, the imaginary time formulation is well-suited and stands as the only simple 
method of calculation. However, for time dependent quantities like two-point correlation functions, the use of this 
formulation requires a continuation to imaginary time and possibly back to real time at the end. On the other hand, 
the real time formulation provides a convenient framework to calculate such quantities, without requiring any such 
continuation at all. 

A difficulty with the real time formulation is, however, that all two-point functions take the form of 2 x 2 matrices. 
But this difficulty is only apparent: Such matrices are always diagonalisable and it is the 11- component of the 
Cji 1 diagonalised matrix that plays the role of the single function in the imaginary time formulation. It is only in the 
calculation of this 11-component to higher order in perturbation that the matrix structure appears in a non-trivial 
way. 

In the literature transport coefficients are evaluated using the imaginary time formulation [6-8] . Such a coefficient 
is defined by the retarded correlation function of the components of the energy-momentum tensor. As the conventional 
perturbation theory applies only to time-ordered correlation functions, it is first necessary to relate the two types of 
OO . correlation functions using the Kallen-Lehmann spectral representation [9-12]. We find this relation directly in real 
time formulation. The time-ordered correlation function is then calculated also in the covariant real time pcrturbative 
framework. 

It suffices to illustrate the procedure with one transport coefficient, say the shear viscosity. It is given by [6, 7, 13], 



V = TO J dhettl J 1 dt ' 1 j d 3 x'9(-t')([7r^(0),7r aP (x',t')}), (1.1) 

where the space integral is over a retarded two-point function. Here (O) for any operator O denotes equilibrium 
. ensemble average, 

Si (0)=Tt(e- pB 0)/Z, Z = Tre- pH , (1.2) 

and Tr a P(x) is the traceless part of the energy- momentum tensor. At temperatures sufficiently below phase transition, 
the pionic degrees of freedom dominate the hadron gas [14]; so one takes only their contribution to this tensor, getting 

n a p(x) = {A p a A' T p -^A a pA pa )d p ${x)-d r7 ${x), A a0 = g af} - u a up , (1.3) 

where <p(x) denotes the pion triplet and is the four-velocity of the fluid, which is (1,0) in the comoving frame. 

In Sec. 2 we derive the spectral representations for the retarded and time-ordered correlation functions in the real 
time version of thermal field theory. The time-ordered function is then calculated to lowest order with complete 
propagators in Sec. 3. We conclude in Sec. 4. 
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FIG. 1: The contour C in the complex time plane used here for the real time formulation. 

II. REAL-TIME FORMULATION 

Here we review the real time formulation of thermal field theory leading to the spectral representations of bosonic 
two-point functions [12]. This formulation begins with a comparison between the time evolution operator e - 8ff (*2-*i) 
of quantum theory and the Boltzmann weight e~^ H = e - lH ( T ~ t t i ~ T ) of statistical physics, where we introduce r as a 
complex variable. Thus while for the time evolution operator, the times ti and t 2 (t 2 > ^i) are any two points on 
the real line, the Boltzmann weight involves a path from r to t — i/3 in the complex time plane. Setting this r = — T, 
where T is real, positive and large, we can get the contour C shown in Fig. 1, lying within the region of analyticity 
in this plane and accomodating real time correlation functions [2, 4]. 

Let a general bosonic interacting field in the Heisenberg representation be denoted by $/(a;), whose subscript I 
collects the index (or indices) denoting the field component and derivatives acting on it. Although we shall call its 
two-point function as propagator, &i(x) can be an elementary field or a composite local operator. (If $;(x) denotes 
the pion field, it will, of course, not have any index). 

The thermal expectation value of the product $;(x)$;,(a;') may be expressed as 

m,n 

where we have two sums, one to evaluate the trace in eq.(1.2) and the other to separate the field operators. They 
run over a complete set of states, which we choose as eigenstates \m) of four-momentum P^, Using translational 
invariance of the field operator, 

$,(a?) = e lP x <P l (0)e- lP x , (2.2) 

we get 

($,(a04(aO) = ^J2 e ~ 0Em ei(km ~ knH *~ X Vl^WlnXnl^]", (0)|m) . (2.3) 

rn,n 

Its spatial Fourier transform is 

d»a; C -*-(«»-«')($ I (x)$j',(a/)} 

^J2 e ~ PEm e l(B ™- B " )(T - r V(fc m - K + fc)H$,(0)|n>H4( )H . ( 2 - 4 ) 



ni.n 



where the times r, r' are on the contour C. We now insert unity on the left of eq. (2.4) in the form 

/CO 
dk' 5(E m -E n + k' ) . 
-oo 

(We reserve ko for the variable conjugate to the real time.) Then it may be written as 

J Sxe-^-^^tix^lix')) = J ^e-<^- T,) M+{k^k), (2.5) 
where the spectral function M + is given by \k' = (k' Q , k)] 

M+(k') = VJe"^ 6 4 (k m - k n + fc')(m|$i(0)|n)(n|4(0)|m) . (2.6) 
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In just the same way, we can work out the Fourier transform of (<&j,(x')<&i(x)) 



tP xe -*»(«-»')(& lt ( a J)$ l ( x ))= I ^O e -<(r-r ') M ~(k' ,k), (2.7) 



with a second spectral function M is given by 



M a ,(k') = ^lj2 e ~ ffEm ~km + k')(m\4(0)\n)(n\^(0)\m) . (2.8) 

The two spectral functions are related by the KMS relation [15, 16] 

M+(k)=e^M^,(k), (2.9) 

in momentum space, which may be obtained simply by interchanging the dummy indices m, n in one of M± (k) and 
using the energy conserving ^-function. 

We next introduce the difference of the two spectral functions, 

Pw (k)^M+(k)~M u ,(k), (2.10) 

and solve this identity and the KMS relation (2.9) for Mj^,(k), 

M+{k) = {1 + f(k )} PH ,(k) , M w (k) = f(k )p H ,(k) , (2.11) 

where f(ko) is the distribution-like function 

■f( fc °) = & 8k!_ i > -oo<fc <oo. (2.12) 

In terms of the true distribution function 

<N) = i^up (2 - 13) 

it may be expressed as 

f(k ) = f(k a ){9(k a ) + 6(-k a )} 

= ne(k ) - 6(-ko) . (2.14) 

With the above ingredients, we can build the spectral representations for the two types of thermal propagators. 
First consider the time- ordered one, 

-iD„,(x,x') = {T * l (x)*l(x')) 

= 6 c (t-t')($i{x)$Ux')) + 6 c (t' -t)($1(x')$i(x)). (2.15) 

Using cqs. (2.5, 2.7, 2.11), we see that its spatial Fourier transform is given by [2] 

( Ll pw {k' ,k)e-^-^{9 c {T - t 1 ) + f(k' )} , (2.16) 

-oo 2lT 

As T — > oo, the contour of Fig. 1 simplifies, reducing essentially to two parallel lines, one the real axis and the other 
shifted by — z/3/2, points on which will be denoted respectively by subscripts 1 and 2, so that t\ = t, T2 = t — i/3/2 
[4]. The propagator then consists of four pieces, which may be put in the form of a 2 x 2 matrix. The contour 
ordered 9's may now be converted to the usual time ordered ones. If t, t' are both on line 1 (the real axis), the 
t and t ordcrings coincide, 6 c (t\ — t{) = 8(t — t'). If they are on two different lines, the r ordering is definite, 
6 c (ti— T2) = 0, 6 c (t2~- t[) = 1. Finally if they are both on line 2, the two orderings are opposite, 9 c (t2~ r^) = 8(t' — t). 

Back to real time, we can work out the usual temporal Fourier transform of the components of the matrix to get 

D u ,(k ,k)= -ap u ,(A£,fc)A(*6,*d), (2.17) 
J-00 27T 
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where the elements of the matrix A are given by [12] 

A n = -A* 2 = 1 + 2mf(k' )5(k' - k ) , 

A12 = A a i = 27Tie l3k 'o/ 2 f(k' )S(k' Q - k ) . (2.18) 
Using relation (2.14), we may rewrite (2.18) in terms of n, 

Al1 = = ¥ k = lV e(k ) + - > 



Ai 2 = Aai = 2niy/n{l + n)e{ko)6(k' - k ) . (2.19) 
The matrix A and hence the propagator Dw can be diagonalised to give 

D w (k ,k)=u( °"' ^*)u, (2.20) 



~D W 



where Du> and U are given by 



D w (k ,k)=l ^M- , f = (^ - ). 

, 2ir k' Q - k - irje(ko) \ \fn VI 



n 



Eq. (2.20) shows that D can be obtained from any of the elements of the matrix D, say D\\. Omitting the indices 
W, we get 

ReD = ReDn , IniD = tanh(/3|fc |/2)ImDn . (2.22) 

Looking back at the spectral functions Mjt defined by (2.6, 2,8), we can express them as usual four-dimensional 
Fourier transforms of ensemble average of the operator products, so that pu> is the Fourier transform of that of the 
commutator, 

pu>(ko,k)= /d*»e tt ^-^<[#,(y),*J'(»')]>. (2-23) 



where the time components of y and y' are on the real axis in the r-plane. Taking the spectral function for the free 
scalar field, 

pa = 2ne{k )S(k 2 - m 2 ) , (2.24) 

we see that D becomes the free propagator, D(ko, k) — — l/(fc 2 — m 2 ). 
We next consider the retarded thermal propagator 

D*( X .x')=i6 c (T-T')({$ l (x,T),$ l ,(x',T')]), (2.25) 

where again r , r' are on the contour C (Fig. 1). Noting eqs. (2.5, 2.7, 2.10) the three dimensional Fourier transform 
may immediately be written as 

Df v (r-r',k) = iB c {t - r') J ^ e ~^-r ) pw ^ fe) _ (2 26) 

As before we isolate the different components with real times and take the Fourier transform with respect to real 
time. Thus for the 11-component we simply have 

r°° dk' 

D*{t-t',k) ll =i0{t-t!) ^e- ik ^-^p u ,{k' Q ,k), (2.27) 



whose temporal Fourier transform gives 



D u ,(k ,k) n - I (2.28) 



o 



l-q 

FIG. 2: The first term in the so-called skeleton expansion of the two-point function. Heavy lines denote full propagators. 
This 11-componcnt suffices for us, but we also display the complete matrix, 

D w (k ,k) = ^ pw{k){ ^: e{ko) + ^e(-k )} -D*?(k) n ) ■ (2 - 29) 

Though we deal with matrices in real time formulation, it is the 11-component that is physical. Eqs. (2.21) and (2.28) 
then show that we can continue the time-ordered two-point function into the retarded one by simply changing the ie 
prescription, 

D$ (fc + «?, fc)n = Du> (ko + ive(qo) -> fc + i% k) . (2.30) 

The point to note here is that for the time-orderd propagator, it is the diagonalised matrix and not the matrix itself, 
whose 11-component can be continued in a simple way. 

III. PERTURB ATI VE EVALUATION 

Clearly the spectral forms and their inter-relations derived above hold also for the two-point function appearing in 
eq. (1.1) for the shear viscosity. We begin with four-dimensional Fourier transforms. To calculate the 11-clcment of 
the the retarded two-point function 

nftfa) =i J d 4 xe^ x -^9(t~t')([7r aP (x,t),n a ^x',t')]), (3.1) 

we consider the corresponding time-ordered one, 

U n (q) = i [ d 4 xe^ x - x '\TTr aP (x,t),7r a P(x' ,t')) , (3.2) 



which can be calculated perturbatively. To leading order, it is given by Wick contractions of pion fields in Tr a p given 
by eq. (1.3). In the so-called skeleton expansion, these contractions are expressed in terms of complete propagators 
(see Fig. 2) to get, 

nii(«) = i J j0^N(l,q)D 11 (l)D 11 (l - q) , (3.3) 
where N(l,q) is determined by the derivatives acting on the pion fields, 

N(l, q) = -6 f(l - qf + ±{l ■ (I - q)} 2 }. (3.4) 

To work out the Iq integral in eq. (3.3), it is more convenient to use An given by eq. (2.18) than by eq. (2.19). Closing 
the contour in the upper or lower half /o-plane we get 

/d 3 l r dW dh" 

J ^p(K,i)^rp(K,i- q )K( qo ,k' ,K), (3-5) 



where 



K _ {i + f(K)}f(K) /(fcoHi + /(*#)} . (36) 

qo - (K - k'o) +ir) q - (k' - k'o) - if)' 
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The imaginary part of Tin arises from the factor K, 

ImK = -^[{1 + /(^)}/(^0 + /(^){1 + /(^0}]%o-(^-^)) 
= -7rcoth03g /2){/(A#) - f(k' )}S(q - (fc - k»)) , 



(3.7) 



while its real part is given by the principal value integrals. 

Having obtained the real and imaginary parts of Hn(q), we use relations similar to eq. (2.22) to build the 11-element 
of the diagonalised IT matrix, 



n = 



d 3 l 
(2tt) s 



dk, 



N(l,q) -±p(k' ,l) -±p(kS,l-q 



2- 



,{1 + /(*&)}/(*&') -/(*o){l + /(*&')} 



2tt 



9o - [k' - kg) + ir]e(q ) 



Finally II can be continued to II ^ by a relation similar to eq. (2.30), 



11 n 



d 3 l 



dk' 



dk'n 



N(l,q) / -p.p%,l)^P-p(kS,l-q) 



2n 



2tt 



,{l + /(^)}/(fcg)-/(^){l + /(fcg)} 
9o - (k' - fc ' ) + iv 



(3.8) 



(3.9) 



Note that in eqs. (3.8,3.9) we retain the /(&o)/(^o) t erms m the numerator to put it in a more convenient form. 
Change the signs of k' and fc ' in the first and second term respectively. Noting relations like 1 + /(— ko) = —f(ko) 
and p(— kg) = —p(ko) we get 



d 3 i 



:N(l,q) 



~2tt^2tt 



p(k ,l)p(k' \l- q)f(&)f(k%)W(q , k' + fc ') , 



where 



W 



1 



qo + k' Q + k'o +irj q - {k' Q + fc ') + it] 



(3.10) 



(3.11) 



Returning to the expression (1.1) for 77, we now get the three-dimensional spatial integral of the retarded correlation 
function by setting q = in eq. (3.1) and Fourier inverting with respect to qo, 

1 J d*x' 6{-t'){[^{Q,Q),n a p{ X ' ,t')\) = - J dq e i «° t 'llK( qo ,q = 0). (3.12) 

This completes our use of the real time formulation to get the required result. The integrals appearing in the expression 
for 77 have been evaluated in Refs. [6, 7], which we describe below for completeness. 

As shown in Ref. [6], the integral over t\, t' and qo in eqs. (1.1) and (3.12) may be carried out trivially to give 



(3.13) 



90=0 



The go dependence of n{\ is contained entirely in W, 



dq 



W(q ) 



90=0 



(k' + k' ' - ir,)* {k' + K + i V f 



2TTid'(k' + fc Q ') . 



(3.14) 



Changing the integration variables in eq. (3.10) from k' , k' ' to fco = k' + k ' and fco = \ {k' Q — fc ') we get 



dH 
(2tt)* 



N(l) 



dko 



where 



(3.15) 



F(ko,l) = 



dko 



ko 



k =0 



(3.16) 



It turns out that the integral over fco becomes undefined, if we try to evaluate F(ko) with the free spectral function 
Po{k) given by eq. (2.24). As pointed out in Ref. [6], we have to take the spectral function for the complete propagator 
that includes the self-energy of the pion, leading to its finite width T in the medium, 



p(k ,l) = - 
1 



1 



1 



(fco - ir) 2 - a; 2 (fc + ir) 2 - uj 2 



= Vl 2 + m 2 



(3.17) 
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Then F(ko,l) becomes 

J.2 J3k o F 2 
( e /3feo_l)2{( fco _ ir )2_ w 2}2{( fco+ir )2_ w 2}2 ' ^ V 

having double poles at fco = 2mn/ j3 for n = ±1, ±2, • • • and also at fco = ±w ± «r. The integral over fco may now be 
evaluated by closing the contour in the upper/lower half-plane to get 

-f- oo j i -I a 

-^F(k , 1) = P n (u,){l + n(u>)} , (3.19) 

where we retain only the leading (singular) term for small T. In this approximation eq. (3.15) gives [7] 

The width T(l) at different temperatures is known [17] from chiral perturbation theory [18], using which r\ has been 
evaluated numerically [7]. 



IV. CONCLUSION 



Here we calculate a transport coefficient in the real time version of thermal field theory. It is simpler to the imaginary 
version in that we do not have to continue to imaginary time at any stage of the calculation. As an element in the 
theory of linear response, a transport coefficient is defined in terms of a retarded thermal two-point function of the 
components of the energy-momentum tensor. We derive Kallen-Lehmann representation for any (bosonic) two-point 
function of both time-ordered and retarded types to get the relation between them. Once this relation is obtained, 
we can calculate the retarded function in the Feynman-Dyson framework of the perturbation theory. 

Clearly the method is not restricted to transport coefficients. Any linear response leads to a retarded two-point 
function, which can be calculated in this way. Also quadratic response formulae have been derived in the real time 
formulation [19]. 
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